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We investigate the entanglement properties of multi-mode Gaussian states, which have some 
symmetry with respect to the ordering of the modes. We show how the symmetry constraints 
the entanglement between two modes of the system. In particular, we determine the maximal 
entanglement of formation that can be achieved in symmetric graphs like chains, 2d and 3d lattices, 
mean field models and the platonic solids. The maximal entanglement is always attained for the 
ground state of a particular quadratic Hamiltonian. The latter thus yields the maximal entanglement 
among all quadratic Hamiltonians having the considered symmetry. 



Classically as well as quantum mechanically the global 
ordering or symmetry of a system often imposes highly 
non-trivial constraints on its local properties. These 
kinds of frustration effects lie at the heart of ordered 
interacting systems, and physicists are faced with these 
phenomena whenever dealing with lattice systems or 
molecular structures. The present paper is devoted to 
investigate how the entanglement of two subsystems of a 
larger system is constraint by such a global symmetry for 
some particularly interesting class of states, the so-called 
Gaussian states Q. 

Gaussian states appear very naturally in several 
branches of physics where entanglement plays a predom- 
inant role. The electromagnetic field in most quantum 
optical setups, atomic ensembles interacting with such 
fields the motion of a collection of trapped ions, or 
the low energy (bosonic) excitations of many interacting 
systems can be very well described by these states. This 
is due to the fact that quantum field theories can be, in 
some regimes, approximated by Hamiltonians which are 
quadratic in some bosonic operators, and thus in thermal 
equilibrium as well as a result of the dynamics the corre- 
sponding states are Gaussian. Thus, there is a growing 
interest in understanding the entanglement properties of 
these states HES0I10 

Our results quantify a very intuitive property of en- 
tanglement, which distinguishes it from the usual corre- 
lations found in classical systems: one particle can share 
entanglement only with a limited number of other par- 
ticles which in turn becomes smaller and smaller as 
the amount of entanglement increases. Furthermore, the 
entanglement that can be shared by a subset of parti- 
cles strongly depends on the symmetries of the multi- 
particle state. For example, if we have a set of particles 
distributed on a lattice in a state with translational sym- 
metry, the maximal entanglement between any two near- 
est neighboring particles should depend on the number 
of spatial dimensions, and should decrease if the total 
number of particles increases. For one-dimensional rings 
of spin i particles quantitative investigations of this kind 
were started in |£g. However, the involved optimization 
problems are highly non-trivial such that up to now only 
a lower bound for the achievable Entanglement of Forma- 



tion (EoF) |lfj is known. In the case of Gaussian states, 
the situation can become even more intriguing since for 
two modes only, the amount of entanglement becomes 
unbounded. If we consider three modes, and impose that 
the global state is invariant under permutations, it turns 
out that the maximum EoF between any pair of modes 
becomes finite. By increasing the number of modes and 
imposing different symmetries to the global state, this 
quantity experiences strong modifications. In this work 
we determine the Gaussian state of N modes which gives 
rise to the maximal EoF between a selected pair of modes, 
for any N and a large variety of symmetry groups. 

From our analysis it also follows that the state for 
which the maximum entanglement is generated under a 
given symmetry corresponds to the ground state of a par- 
ticular Hamiltonian, quadratic in the bosonic operators, 
which can be easily constructed. Thus we can deter- 
mine the Hamiltonian, invariant under a certain symme- 
try group, that generates the maximum two-mode en- 
tanglement for the physical systems mentioned above. 

Although we will consider rather general symmetry 
groups, we will illustrate our results for groups which 
can be associated to symmetric graphs (Fig^a)) , since 
they give an intuitive geometric depiction of the group 
and they are the ones that naturally appear in many 
physical systems. For example, we will give the optimal 
EoF for states that have the symmetries of a lattice in 
any dimension, including square, cubic, hexagonal, and 
trigonal lattices (Tab 01, or those of all platonic solids 
(Tab©. 

Let (Qi, . . . , Qn,Pi, • ■■> Pn) ■— R be the N conjugate 
pairs of canonical operators characterizing N modes and 
obeying the canonical commutation relations [Rk,Ri] = 
icjfcz with a being the symplectic matrix (cf.[l|). Let us 
consider a subgroup G of the permutation group and two 
particular modes, k, I < N, for which there exists a group 
element such that g(k) = I and g(l) = k (this condition 
will be relaxed later on). We construct a Hamiltonian 
operator as follows 

Hmax = TjTji (Qk*0 + Qg(')) 2 + ( p g(k) - Pgii))^ 1 ) 
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Let us denote by Eq the ground state energy and by 
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FIG. 1: (a): Apart from chains, cubic lattices and meanfield 
clusters there are several familiar symmetric graphs. Exam- 
ples are the five platonic solids (e.g. the dodecahedron), and 
hexagonal or trigonal lattices, (b): Maximal nearest neighbor 
entanglement E max (ebits) in a ring of N harmonic oscillators. 
The dotted curves represent the envelopes corresponding to 
Ea.ltTtjl. 



'Lo the corresponding ground state, which is a Gaussian 
state, i.e., has a Gaussian Wigner function. The rela- 
tion between this Hamiltonian and the EoF will later on 
be established by linearizing the expression for the latter. 
We will show that the Gaussian state which is invariant 
under G and which maximizes the EoF of the modes k 
and I is exactly an d that the corresponding EoF is 
E m ax = E F (E ), where 



E F (A) ■- 
and c±(A) 



C+ (A)log[c + (A)]- C _(A)log[c_(A)], (2) 



(A -1 / 2 ± A 1 / 2 ) 2 /4. Hence, the ground 
state of H max bas maximum entanglement under all 
G-invariant Gaussian states and conversely, among all 
quadratic G-invariant Hamiltonians H max generates the 
largest amount of entanglement at zero temperature. 

The ground state of Hamiltonians of the form Q can 
be easily determined by resorting to symplectic space. 
We define the covariance matrix (CM) of a Gaussian state 
p as usual, 



T kl := ({{R k - (R k )), (R k - (R k ))}+), 



(3) 



which must fulfill r > ia Q. Let us also introduce the 
Hamiltonian matrix |f| corresponding to H max as H = 
H+ © H_ with [HI 
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(4) 



tk,l) 



±[\k)(k\ + \l)(l\±(\k)(l\ + \l)(k\)]. (5) 



The matrix H > can be diagonalized by a symplectic 
matrix S, H = SDS T . Since tr \pH~\ = tx\TH], the 



lattice 


Emax 


N a 


hexagonal (2d) 


10.61 


3 


square (2d) 


6.31 


4 


trigonal (2d) 


2.69 


6 


cubic (3d) 


2.62 
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TABLE I: Maximal nearest neighbor entanglement E max (in 
units of 10 -2 ebits) for some infinite 2d and 3d lattices. N a 
is the number of adjacent vertices. 

ground state energy of H is given by [l4j 
Eq 



inftrpTffl 
r 



inftr [YD] 



h] /2 h_h] /2 \ 



(6) 

Since H has a null space, the CM Tq of the ground state 
>Fo of H is achieved in the limit e — > of T e := S^^S^ 1 , 
where S e is the symplectic matrix that diagonalizes H + 
el. 

Before we prove the above statements, let us utilize the 
results to analyze the maximum EoF for several interest- 
ing symmetry groups. We will concentrate on groups 
which can be associated to a symmetric graph. 

Consider a simple undirected graph with N vertices, 
characterized by an adjacency matrix A, which is such 
that A k i = 1 if the vertices k and I are connected by an 
edge, and A k i = if there is either no edge or k = I. The 
symmetry group G of the graph contains all permutations 
g, which commute with the adjacency matrix [A, g] = 0. 
The graph is called symmetric if all edges as well as all 
vertices are equal in the sense that every edge and every 
vertex can be mapped onto every other one by an element 
of G. Examples of symmetric graphs are given in Fig^a) 
and Tab lfllll By utilizing this symmetry we can simplify 
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(7) 



where £ = {(k,l)\A k i = 1} is the set of edges. Hence, 
the sum in Eq.Q) runs over all edges which correspond 
therefore to physical interactions between the adjacent 
modes. By observing that [H + ,H_] = we can derive 
the ground state CM (or the respective limit): 



r = \IH-H. 



h+h: 



(8) 



Note that when acting on two modes only, the ground 
state of H max is the original singular EPR-state. More- 
over, since i7_ has a Kernel containing the vector 
(1, 1, . . . , 1) the ground state of any H max will always 
be singular and the maximal entanglement is thus only 
attained exactly in the limit of infinite squeezing. 

In the following we will apply the obtained results to 
some examples of familiar symmetric graphs. 

Chains and rings: The simplest non-trivial exam- 
ple of a symmetric graph is a ring of N nodes represent- 
ing translation and reflection symmetry. In this case the 
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operators H± have the form 
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H± = —[21±(T + T- 1 )] 



(9) 



where Tjy = 6 k ; + i mo djv ^ s the cyclic shift operator. H± 
are so-called circulant matrices which can be diag- 
onalizcd simultaneously by a Fourier transform yielding 



1=0 



N 



cot 



AT' 



TV even 



jj cot ^ , N odd 



(10) 



Hence, the entanglement remains finite and is suppressed 
in rings with an odd number of modes (see Fig^l. It 
approaches E max = 0.30 ebits (Eq = -) in the limit 
N — > oo. This value is comparable to the 0.29 ebits 
conjectured for an infinite chain of spin | particles 0- 

Cubic lattices: The result obtained for the chain 
has a straight forward extension to <i-dimensional cubic 
lattices. Imposing periodic boundary conditions (i.e. a 
lattice on a torus) for a cubic lattice of N d modes we get 



4N d 



21 ± 



o=l 



(11) 



where now T( Q ) is the shift operator acting on the a'th 
of d tensor factors, each corresponding to one of the di- 
mensions of the lattice (e.g. T( 2 ) = 1 <8> T ® 1 . . .). Diag- 
onalizing H± by a tensor product of Fourier transforms 
leads to 



En — 



N N 
J\fd ' ' ' 



h=l 



ld=i L 



1 2tt x- 



fc=l 



, (12) 



which goes to 1 (E max — > 0) for e? — > oo, and is calcu- 
lated explicitly in Tab[I] for the infinite two and three 
dimensional cubic lattice. 

Mean field clusters: When every mode is con- 
nected to every other one, i.e., when we have complete 
permutation symmetry, then 



H± = [2N(N - l)]" 1 [(AT - 1)1 ± (E — 1)] . 



(13) 



where E^ = 1, which leads to Eq = y ^j^-- Hence, the 
maximal entanglement decreases with the number N of 
modes and vanishes as ~ log N in the limit — * oo. 

Platonic solids: The results for the graphs corre- 
sponding to the three dimensional platonic solids can be 
found in Tab.(|n|). 

All these examples indicate three different tendencies 
for the maximal EoF: 

1- E max decreases with the number of adjacent ver- 
tices. 



platonic solid 


Emax 


N a 


JV 


Eo 


tetrahedron 


19.74 


3 


4 


1 

72 


cube 


19.74 


3 
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1 

V2 


dodecahedron 


11.12 


3 


20 


^(12 + 5^ + 2^) 


octahedron 


10.75 


4 


6 


|(3 + V3) 


icosahedron 


5.37 


5 


12 


— + 4= 



TABLE II: Maximal amount of nearest neighbor entangle- 
ment Emax (measured in units of 10 -2 ebits) and the respec- 
tive ground state energy (minimal EPR-uncertainty) _Eo for 
the five platonic solids. N a is the number of adjacent vertices 
and iV the total number of nodes. 

2. E max decreases with the total number of vertices. 

3. E max is suppressed in loops with an odd number of 
vertices, which give rise to additional frustration. 

Let us now to proceed to prove our main result. We 
denote by F and 7 CM of the global state and the reduced 
density operator for the modes fc, I'm whose entanglement 
we are interested. The first CM must fulfill 
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,T g )T(T g ®T g ) T . 



(14) 



The CM 7 of a two-mode subsystem can always be 
written, up to local symplectic transformations Si$, in 
the standard form Q 



(Si © S 2 ) 7 (5i © S 2 




(15) 



The fact that G contains by assumption an element which 
maps k <-> I immediately implies that tia — tlb and Si = 
S2 = S. Hence, given a global CM we can always find 
another one given by ( 0^Li S)r( i=1 S) , which is 
also G-symmetric, and such that 7 has the standard form 
(I15JI with ua = Tib =:n. 

We are interested in maximizing the EoF of 7. Since 
riA = we can use the results of 0, which show that 
this quantity is given by Ep(A), where the function Ep 
has been given in J2J, which is a monotonically decreas- 
ing function of the so-called EPR-uncertainty A. Thus, 
maximizing the EoF is equivalent to minimizing A. This 
last quantity is a highly nonlinear function of the param- 
eters n, k qi k p , and thus minimizing it with respect to all 
possible global V looks as a very daunting task. In order 
to overcome this problem, the trick is to linearize the ex- 
pression of A by including an extra maximization in the 
problem, i.e. writing 



A = inftr 

s>0 



inf tr 

s>0 



l h (k,D 
s 



7 [shy 

V s 



(16a) 
(16b) 



In the last step we have used that 7 is the reduced CM 
of T, and Ea. (|Tljl . 
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We show now that T must correspond to a pure state. 
If a CM r m corresponds to a G-symmetric mixed state, it 
can always be decomposed into a G-symmetric pure state 
CM T„ and a G-symmetric matrix M > via T, n — T p + 
M [l5| . This decomposition can be interpreted as adding 
classical Gaussian noise to T p £| . Since this will certainly 
not increase the entanglement between any two modes of 
the system, maximal entanglement will be attained if the 
overall CM T corresponds to a pure state. 

We can thus exploit the fact that every pure state CM 
can be written as 

r = ( X XY , ) , (17) 

^ YX YXY + X- 1 J 

with X > and Y = Y T 0. Hence, the EPR- 
uncertainty l|lrja|) becomes 

A = inf str [XH+] + -tr [(X" 1 + YXY)HJ\ , (18) 

Maximizing the entanglement means minimizing A with 
respect to X and Y under the constraint that they pa- 
rameterize a G-symmetric CM. We can, however, drop 
this constraint since the symmetry of H± will force the 
optimal To to have the right symmetry. Moreover, Tq 
will be the ground state corresponding to the Hamilto- 
nian matrix H — H + © H— since 

A = inf A = inf tr [XH+] + tr [X^HJ] (19) 
= inftr[r(ff+ ®H-)\ (20) 

where we have first set Y = 0, since tr [YXYH-] > 
and then incorporated the infimum over s into that over 
X. This completes the proof of our main result, since Ao 
is by Eas. (|4l2U[) equal to the ground state energy Eq of 
the Hamiltonian in Eq.Q. 

By imposing less restrictions on the symmetry group 
than requiring the existence of an element which inter- 
changes k <-> Z, the maximal achievable EoF could grow. 
In the examples on chains and lattices, we implicitly im- 
posed the translational and reflection symmetry. The 
same results can however be derived with only the trans- 
lational symmetry. In fact, the presented proof can be 
extended in a straight forward manner to all symmetry 
groups with Abelian commutant, including those of rings 
and cubic lattices, without imposing reflection symmetry. 
The proof can be found in the appendix. 

In conclusion, we have determined the maximal entan- 
glement between two modes under the constraint that 
the overall system is in a Gaussian state which has some 
symmetry with respect to the ordering of the modes. The 
result was derived by linearizing the entanglement func- 
tional which permits to perform the maximization in a 
simple way. We find that the maximal entanglement is 
connected to the ground state of a particular quadratic 



Hamiltonian which possesses the same symmetry as the 
state. The state that maximizes the EoF is precisely the 
ground state of such a Hamiltonian. The maximal entan- 
glement turned out to be finite in all the discussed cases, 
and is even comparable to the (conjectured) values for 
spin I systems for the case of rings 0. We have shown 
how the entanglement decreases with the number of spa- 
tial dimensions, and how it depends on the geometry of 
the state. Finally, although we have concentrated here 
on Gaussian states, we think that similar results may be 
obtained for the case of qubits. In that case it is also 
possible to linearize the expression for the EoF ^| an d 
in this way to relate it to the ground state energy of a 
specific nearest-neighbor interaction Hamiltonian, which 
allowed us to solve the problem for Gaussian states. 
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APPENDIX: GROUPS WITH ABELIAN 
COMMUTANT 

In this appendix, it will be shown that the reflection 
symmetry need not be imposed if the commutant of the 
group is Abelian, as the states with maximal EoF au- 
tomatically have this property. The proof is as follows. 
The reduced covariance matrix corresponding to modes 
k and I of a state with a symmetry that maps g{k) = I is 
of the form 

(a b 

b c 



7 : 



d e\ 

f 9 
d f a b 
\e g b cj 



It can readily be checked that this can be brought into 
symmetric normal form 7 by local symplectic transfor- 
mations of the form 

7 = (o e i)(s © 5)7(5 © sf(o T © 1) 

with 5 symplectic and O orthogonal. Suppose now that 
7 and the global covariance matrix T are such that 5 = 
1, which can always be done by a transformation T 1— » 
(©iI 1 5')r(ffi^ 1 5') T . Then the EPR-uncertainty A is 
given by 



A = inf Tr 

s>0 



7 (O t © 1) 



( s s \ 

l/s -l/s 

s s 

\0 -l/s l/s j 



(0©1) 



Parameterizing O = 



cos(0) -sin(0) 
sin(0) cos(0) 



one gets 



'(/</ 



h — h T 



'I'll 



scos 2 (0) + 5!h!W SC os(0) 
scos(0) s 

(±-s)cos(0)sin(0) 
— s sin(0) 



• 2//i\ 1 cc 

ssm (9) -\ — 

cos(g) 



-(e) 








Let's next extend these Hamiltonian blocks to N- modes 
by averaging over the group: 



Ha/3 = ^2 T 9 ha P T i 

1 1 gee 



The goal is now to calculate the sum of the symplectic 
eigenvalues of this Hamiltonian. As the commutant of the 



group was assumed to be Abelian, all N x TV matrices 



A; 



1 

|G| 



— Vr 9 |z)(j|T ( 
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S {0,1}, commute and hence can be diagonalized 
simultaneously by a unitary operator (eventually com- 
plex), yielding diagonal elements , where G {0, 1} 
and 1 < ji < N . Applying the same unitary transforma- 
tion to the four matrices H a p, the complete Hamiltonian 
becomes a direct sum of 2 x 2 blocks of the form 



B,, 



b^ 

b^ * 

qp 



"qp 

pp 



(24) 



In general, the symplectic eigenvalues of a positive op- 
erator r are given by the square roots of the eigenvalues 
of the operator (roTV T ), and the unitary operation un- 
der consideration does not change these eigenvalues. It 
can now readily be checked that the sum of the symplec- 
tic eigenvalues corresponding to the block is given by 
the formula 



v£ + v» = yjb% q bb - (Re(^ p )) 2 
Due to equation (|24(l . it holds that 



Vt-UL - (Re(6^)) 2 = (A>) T QA> 



where 



A p = 



Q = 




sin 2 (0) 



1 

s 

s 



2/001 

0-2 
\l 



The presence of the second term in Q (which is not there 
when reflection symmetry is present) can never lead to an 
improvement of the entanglement, as it can only increase 
the symplectic eigenvalues. Indeed, it holds that 

^00^11 ~ l^oi I 2 ^ 0, 

as it is the determinant of a principal submatrix of a 
positive 2N x 27V matrix P. (P is obtained by first ap- 
plying the group symmetry to the hypothetical blocks 



h qq = |0)(0|;V = h p q = \0)(M\hpp = 1 1> <1 1 (which 



is manifestly positive), and then diagonalizing the four 
blocks by a unitary transformation.) Therefore, the op- 
timal choice for the orthogonal matrix O is to choose it 
equal to the identity, imposing that the optimal solution 
obeys the reflection symmetry. 



